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abstract 

The  eikonal  equation  of  ray  acoustics  is  discussed  for  the  general 
case,  with  the  index  of  refraction  a  function  of  the  three  Space 
coordinates..  Two  examples  illustrating  the  influences  of  a  two- 
dimensional  heat  source  and  a  thermal  mixing  zone  on  the  acoustic 
ray  paths  are  presented.  Numerical  results  show  that  for  long  range 
acoustical  ray  tracing  inhomogeneities  of  the  ocean  in  the  horizontal 
plane  can  cause  refraction  effects,  which  are  not  negligible  in  com¬ 
parison  with  effects  in  the  vertical  direction.  Rays,  which  travel 
in  the  vicinity  of  a  heat  source,  can  be  refracted  as  much  as  6°. 
Thermal  mixing  zones  can  increase  the  range  of  shadow  zones  near  the 
ocean  surface  up  to  50%.  Even  new  shadow  zones  can  be  formed. 
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1.  Introduction 

In  oceanography  ray  acoustics  is  usually  applied  to  media  in  'which 
the  speed  of  sound  changes'  only  in  one  direction,  [l],  [2].  This 
restriction  is  justified  for  most  practical  purposes  because  variations 
of  the  speed  of  sound  are  much  larger  in  the  vertical  direction  of  the 
ocean  than  in  the  horizontal  plane.  Better  knowledge  Of  the  ocean 
currents  in  recent  years  raises  the  question  about  the.  influence 
of  a  general  three-dimensional  gradient  of  the  speed  of  sound  on  long- 
range  acoustical  ray  tracing. 

For  this  reason  the  Naval  Weapons  Laboratory  has  prepared  a  program 
for  the  solution  of  the  eikonal  equation  in  which  the  speed  of  sound 
can  be  a 'general  function  of  the  space  coordinates.  This  program  will 
be  used  in  the  near  future  to  solve  special  practical  problems. 


In  this  study  the  basic  equations  are  discussed,  and  the  coded 
program  is  applied  on  two  models  representing  a  heat  source,  and  a 
thermal  mixing  zone  in  the  ocean.  These  models  are  rather  artificial 
because  the  detailed  features  of  oceanic  inhomogeneities  are  complex, 
and  transient.  They  Serve,  however,  to  demonstrate  the  magnitude  of. 
refraction  one  may  expect -in  the  oceans'  und>jer  such  conditions.  The 
influence  of  small-scale  fluctuations  has  not  been  considered. 


2.  Basic  Assumptions  and  Equations 

All  ocean  currents  are  very  slow  compared  with  the  speed  of  sound 


so  that  their  velocities  can  be  neglected.  They  are,  however,  often 
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connected  with  large  gradients  of  temperature  and  salinity ,  which 
influence  the  speed  of  sound.  The  medium  in  consideration,  therefore, 
can  he  regarded  as  isotropic  hut  inhomogeneous .  Presuming,  furthermore, 
the  validity  of  ray  acoustics  —  that -the  speed  of  sound  does  not  change 
much  in  distances  of  the  order  of  one  wave  length  —  one  can  use  the 
eikonal  equation  (see  e.  g.  Chernov  [3]): 


where  s  -  arc  length  along  the  ray  path  from  a  given  starting  point 
r  “  ?(s)  the  position  vector  of  the  ray 
n  »  index  of  refraction  of  the  medium 
In  a  cartesian  coordinate  system  (x,  y,  z)  equation  (1)  is  equivalent 
to  the  following  three  ordinary  differential  equations  of  second  order: 


where  the  dots  denote  the  derivatives  with  respect  to  s.  The  ray  path 
is  determined  by  the  initial  conditions 
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The  set  of  equations  (2)  through  (4)  is  programmed  for  solution  by  the 
Punge-Kutta  method  on  the  computer  IIM  -  7090. 

By  elimination  of  the  implicit  parameter  s  equations  (2)  through 
(4)  may  be  combined  into  two  equations: 

+  /  +  (6) 

*'  *  U .+  y?  +  *^4*  «'  -4r|S  >-°  <7> 

The  primes  here  indicate  differentiation  with  respect  to  x,>  The. 
initial  conditions  (5)  become 

x*x0:  y0»  z0J  yo  >  zo  •  (8) 

The  index  of  refraction  is  a  function  of  the  temperature  x,  the 
salinity  S,  and  the  pressure  p  of  the  medium: 


n  -  n(T,  S ,  p) 

(9) 

with 

T  *  T(x,  y,  z) 

(10) 

S  =  S(x,  y,  z) 

(ID 

P  “  P(z) 

(12) 

if  the  coordinate  z  is  chosen  normal  to  the  ocean  surface.  For  the 
calculation  of  the  index  of  refraction  (9)  as  a  function  of  T,  S,.and 
p  the  formula  by  Wilson  [4] .has  been  used.  The  functions  (10)  and  (11) 
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depend  on  the  medium  under  consideration.  p(z)  is  given  by  the  hydro¬ 
static  equation.  From  (9)  through  (12)  there  follows 


n  =  n(x,  y,  z) 


(13) 


Equations  (6)  and  (.7)  are.  in  general  too  complicated  to  yield  an 

analytic  solution  if  the  index  of  refraction  is  a  function  of  more  than 

one  variable.  .Only  in  special  cases  are  analytic  solutions  possible. 

Two  such  cases  are  briefly  examined  for  n  =  n(x,  y)  in  the  next  section. 

3.  Discussion  of  the  Eikonal  Equation  for  the  Two-dimensional  Case 

A  medium  will  be  assumed  in  which  the  index  of  refraction  is  only 

a  function  of  the  coordinates  x  and  y,  i.  e.  -^5  *  0.  In  addition  only 

3z 

those  rays  will  be  considered  which  start  out  initially  in  the  (x,  y) 

plane.  Then,  because  of  the  restriction  —  =  0  they  will  remain  in 

dz 

that  plane.  The  differential  equations  (6)  and  (7)  reduce  to 


A/I  ia  _  I 


y  +  (1  +  y  )(nSy 


n  3y' 


(14) 


with  the  initial  conditions 


x  *  Xg : 


y0»  y©  • 


(15) 


In  general,  the  differential  equation  (14)  cannot  be  solved  in  a 
closed  form.  There  exist,  however,  classes  of  analytic  solutions 
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which  will  be  discussed  briefly.  ‘ 

A  first  integral  can  be  obtained  if  the  last  factor  in  equation  (14) 
represents  an  exact  differential  dY  : 


-  Is  dy  -  i  I?  dx  =  dY  (16) 

n  ax  n  qy  • 

According  to  Kamke  [5],  the  first  integral  is 

y'  =  -  tan(Y  +  const.)  (17) 

Equation  (16)  holds  if  the  function  n  satisfies  the  condition 

nV2n  =  jVn  \s  (18) 

or:  n  «  0!  (x  +  iy)02(x  -  iy)  (19) 


Where  only  real  values  of  n  are  of  interest,  fa  and  fa  are  arbitrary 
functions  of  the  Complex  number  x  +  iy  and  its  conjugate  x  -  iy..  In 
many  practical  cases  the  approximation 

1  is  5  (20) 

n  ox  Ox 

may  be  used.  The  condition  (18),  then,  is  replaced  by  the  Laplace 
equation  V®n  =  0  with  the  well-known  general  solution 
n  =  Cft  (x+  iy)  +  %(x  -  iy)  . 
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Another  class  of  analytic  solutions  of  equation  (14)  can  be  found  by 
choosing  V 

£  y'  -  f  - 0  <21> 

The  ray  paths  are  straight  lines  because  y#  =  0.  Besides  the  trivial 
case  n  -  const,  which  is  also  a  solution  of  (18) ,  equation  (21)  will  be 
satisfied  if 

y '  -  K  -  const.  (22) 

3y  dx 

The  solution  of  equation  (22)  is: 

n  “  n(x  +  Ky)  (23) 

e 

Thus,  in  a  medium  described  by  the  above  index  of  refraction,  one  ray 
(with  the  initial  slope  K)  will  be  propagated  as  a  straight  iine. 

4.  Two-dimensional  Heat  Source  in  an  Infinite  Ocean 

The  ray  paths  near  a  heat  source  will  be  calculated  numerically  by 
considering  only  paths  in  the  horizontal  plane  (x,  y) .  In  practice  a 
heat  source  may  be  represented  by  a  melting  iceberg  or  by  a  separated 
vortex  flow  with  a  temperature  gradient,  as  it  occurs  in  and  near  the 
Gulf  Stream.  The  heat  transfer  is  approximately  described  by  the 
transient-state  law 

T  -  T»  _r  ot  Vs  (24) 

’  —  s 
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where  r*  =  x?  +  ya  .  The  heat  source  is  located  at  the  origin  r  =  0 . 
T  is  the  temperature  in  the  plane  (x,  y) ,  T0  the  temperature  at  the 
origin,,  and  T  .  the  temperature  at  infinity1*.  The  coefficient  or  is 
inversely  proportional  to  the  time.  For  ray  tracing  the  dependency 
on  the  time  can  be  neglected  since  the  time-dependent  heat  transfer 
is  small  with  respect  to  the  speed  of  sound.  Using  a  linear  approxi¬ 
mation  for  the  speed  of  sound  as  a  function  of  the  temperature  one 
obtains 


n  =  1  +  ke"^  r  (25) 


where  k  is  a  constant  proportional  to  the  temperature  difference 
between  the  heat  source  and  the  ocean  at  a  large  distance  from  the 
heat  source. 

Equation  (14)  yields  with  equation  (25) : 


// 

y 


+ 


Cl  +  y/3) 


■3gk(y."  siX 


=  0 


(26) 


Substituting 


§  =  yk  xj.  T1  *  jot  y 


(27) 


one  arrives  at  the  normalized  equation 


if  +  -an  *  it")—  )-  0  <28> 

where  the  primes  denote  here  the  derivatives  with  respect  to  5  . 
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Numerical  solutions  by  means  of  the  IBM-7090  computer  are.  shown 
in  Figures  1  through  4  for  the  following  parameters 
kj.  =  -  0.03155;  (T0  =  10  °C,  T«  -  0°C)  ) 


ks  =  -  0.06310; 


(T0  =  20  °C,  T^  =  0°C) 


heat  source 


k3  -  +  0.03155; 
k4  =  +  0.06310; 


(T0  =  10°C,  T^  =  20 °C) 
(T0  =  0°C,  ^  =  20 °C) 


heat  sink 


For  the  calculation  of  k  a.  linear  approximation  of  ‘Wilson's 
formula  [4]  has  been  applied. 

The  curves  in  Figures  1  through  4  can  be  used  to  construct  the 
ray  paths  from  a  point  source  of  sound.  Two  examples  are  shown  in 
Figures  5  and  6,  which  are  drawn  with  the  aid  of  Figures  2  and  4, 
respectively. 

Discussion  of  the  results. 

The  divergence  of  ray  paths  is  a  measure  of  the  intensity  of 
sound  waves.  It  can  be  seen  that  a  heat  source  causes  the  bundles 
to  diverge  and  thus  diminishes  the  intensity  whereas  a  heat  sink 
has  the  opposite  effect.  Besides  this  qualitative  statement  the 
total  deviation  of  the  ray  paths  can  be  found.  Figure  7  shows  the 
total  angle  of  deviation  A0  as  a  function  of  T|(-co)  —  the  perpen¬ 
dicular  distance  between  the  asymptote  of  the  ray  and  the  heat 

source  —  for  T0  -  T_  ■=■  20°C.  The  solid  curve  is  obtained  from 
w  00 
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the  numerical  solutions ,  and  the  dashed  curve  is  the  graph  of  the 
function  • 

A9  =  2Vrr  k  11(-®)e"^^"'8^8  (29) 

Equation  (29)  is  derived  in  Appendix  A  as  a  first  order  approxima¬ 
tion  without  knowing  the  actual  ray  paths.  The  agreement  of  the 
two  curves  is  satisfactory. 

From  Figure  7  it  is  seen  that  rays  passing  near  a  heat  source 
can  be  refracted  up  to  6°.  This  may  be  regarded  as  a  measure  for 
the  magnitude  of  error  that  can  occur  if  a  change  of  the  refraction 
index  in  the  horizontal  direction  of  an  ocean  with  a  heat  source 
is  not  considered. 

5.  Thermal  Mixing  Zone  in  an  Ocean 

In  this  example  ray  tracing  in  the  vertical  cross-section  of 
the  boundary  layer  between  two  fluid  flows  normal  to  the  flow 
directions  is  examined.  The  flows  have  different  temperatures  so 
that  a  thermal  mixing  zone  associated  with  the  momentum  boundary 
layer  exists. 

Let  the  cross-section  be  the  (y,  z)  plane  which  is  normal  to 
the  flow  direction  x.  The  horizontal  temperature  distribution  at 
the  ocean  surface  z  =  0  in  the  y-direction  is  assumed  to  be  propor¬ 


tional  to  the  flow  velocity  and  according  to  Schlichting  [6]  : 
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where  T(-“,  0)  and  T(“,  0)  are  the  surface  temperatures  of  the 
colder  fluid  and  the  warmer  fluid,  respectively,  and  y  is  the 
point  of  inflection  of  T(y,  0).  For  the  vertical  temperature  distri¬ 
bution  of  the  colder  fluid  an  empirical  function  h(z)  is  used,  which 
is  typical  for  the  ocean  (see  Fig.  8) [2]: 


h(z) 


T(-m,  z)  -  T(-«,  0) 
T(-“,  0)  -  T(-“,  zB) 


(31) 


where  T(-°°,  zb)  is  the  temperature  of  the  colder  fluid  at.  the  bottom. 

The  function  h(z)  is  obtained  by  linear  interpolation  of  the 
temperature  vs.  depth  graph  in  Figure  8.  The  points  between  which 
this  interpolation  is  made  are  tabulated  in  Appendix  B.  The 
temperature  is  constant  for  0  s  z  £  60.  The  main  thermocline 
extends  from  z  =  60  to  about  1500  m.  For  depths  greater  than  3000  m, 
the  temperature  is  constant  again.  Figure  8  also  shows  the  speed 
of  sound  V  vs.  depth  for  the  colder  fluid  (y  =  -») .  It  will  be  ob¬ 
served  that  V  has  a  relative  maximum  at  z  =  60  m.  The  minimum 
occurs  at  z  =  600  m,  which  is  the  depth  of  the  deep  sound  channel. 

In  order  to  obtain  the  temperature  distribution  over  the  entire 
(y,  z)  plane  it  is  assumed  that  for  any  line  z  =  const  the  tempera¬ 
ture  distribution  is  similar  to  that  for  z  =  0  and  that  the  dif¬ 
ference  T(-“>,  z)  -  T(+»,  z)  varies  in  the  same  way  as  h(z).  The 
complete  temperature  distribution  in  the  (y,  z)  plane  can  then  be 
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expressed  by 

T(y ,  z)  =  T(-«,  0)  +  0)  -  T(-»,  zB)]h(z) 

+[T(®,  u)  -T(-«,  0)  ][1  -f  h(z) ]erf  (y  -  y)  .  (32) 

Assumption  (32)  implies  that  the  velocity  of  the  flow  as  a  function 
of  the  depth  is  proportional  to  the  vertical  temperature  distribu¬ 
tion.  This  hypothesis,  although  not  borne  out  conclusively,  is  at 
least  not  contradicted  by  available  experimental  results.  See 
Stommel  j*7  ]  . 

The  surface  of  the  ocean  is  taken  as  a  specular  reflector,  and 
the  bottom  is  assumed  to  be  a  perfect  absorber. 

For  numerical  calculations  the  following  data  are  chosen  which 
represent  typical  values  for  the  Gulf  Stream: 

T(-%  0)  =  12°C,  T(+®,  0)  =  27 °C 

T(-°°,  z8)  =  2°C 
Zg  =  5000  m 
S  =  35  °/00 

p  =  (1.033  +  0. 10 23 z) kg  cm”2 

The  computed  ray  paths  are  shown  in  figures  9  through  17. 

Each  ray  ib  identified  by  the  angle  it  makes  initially  with  the 
ocean  surface.  The  sound  source  is  located  at  three  different 
depths:  z  =  0,  600  and  3000  m  for  each  of  two  positions  Inside 
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the  thermal  mixing  zone:  at  y  =  0,  and  y  =  30  km.  The  latter 
position  corresponds  to  the  inflection  point  y  of  the  error  func¬ 
tion  (30).  The  scale  of  y  is  then  determined- by  letting  the 
distance  between  y  =  0  and  y  =  30  km  correspond  to  three  standard 
deviations  of  the  error  function  (30).  Computations  were  carried 
out  up  to  a  range  of  100  km. 

Discussion  of  the  Results 

Figures  9  through  11  show  ray  paths  in  an  ocean  in  which  the 
speed  of  sound  varies  only  with  depth  as  in  Figure  8. 

With  the  sound  source  at  the  surface  (Figure  9)  all  ray?  with 
an  angle  of  inclination  0O  <  2°5/  are  confined  in  a  shallow  layer 
60  m  below  the  surface.  Rays  with  2°5  '  £  0O  ^  13°55  ‘  are  bent 
downwards '  till  they  reach  the  depth  of  the  sound  channel  (600  m)  , 
and  then  are  bent  upwards  again.  This  results  in  a  large  shadow 
zone  extending  up  to  a  range  of  47  km  -  region  A.  Rays  with 
0O  >  13°55  /  are  absorbed  by  the  bottom,  so  that  two  other  shadow 
zones  —  regions  B  and  C  —  are  formed.  If  the  bottom  were  removed 
(ze  =  <»)  these  shadow  zones  would  not  exist,  but  one  can  easily  see 
that  the  intensity  would  be  relatively  low. 

In  Figure  10  the  sound  source  is  placed  at  the  deep  sound  \ 
channel  depth  (z  =  600  m) .  The  rays  oscillate  about  this  channel, 
which  is  a  well  known  result  [2],  [8]  .  Three  shadow  zones  are 
formed:  regions  A,  B,  and  C.  In  contrast  to  Figure  9  the  regions 


12 


_ _ _ _ _  NWL  REPORT  NO.  1854 

A  and  B  extend  up  to  the  surface. 

In  Figure  11  the  sound  source  is  located  3000  m  below  the  surface, 
at  the  depth  where  the  vertical  change  of  the  temperature  becomes 
zero.  All  rays  with  60  s  12°  are  bent  up  to  and  reflected  by  the 
surface,  those  with  90  >  12°  are  absorbed  by  the  bottom,  causing  the 
formation  of  two  shadow  zones:  regions  A  and  B. 

Figures  12  through  14  show  the  ray  paths  in  an  ocean  with  the 
two-dimensional  temperature  distribution  (32) ,  and  are  directly  com¬ 
parable  with  the  preceding  Figures  9  through  11.  In  general  it  can 
be  noted  that  the  effect  of  the  thermal  mixing  zone  is  to  depress 
the  rays  downwards ,  and  thus  to  enhance  all  surface  shadow  zones , 
and  even  create  some  new  ones. 

A  comparison  of  Figures  9  and  12  shows  that  the  shadow  zone  A 
now  extends  up  to  a  range  of  60  km  instead  of  47  km  —  an  increase 
of  29% .f 

The  influence  of  the  thermal  mixing  zone  has  created  another 
interesting  effect,  which  is  not  shorn  in  Figure  12.  Because  the 
rays  are  gradually  shifted  downwards,  all  of  the  rays  with 
90  <  2°5  '  axe  not  confined  any  more  within  the  shallow  60  m  surface 
layer  as  in  Figure  9.  Every  time  the  limiting  ray  in  this  surface, 
layer  reaches  its  n'&dir,  a  small  packet  of  rays  will  have  been 
shifted  downwards  far  enough  so  that  it  penetrates  the  60  m  depth 
and  is  subjected  to  the  negative  speed  of  sound  gradient  of  the 
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main  thermocline.  To  show  this  effect  the  ray  paths  for  five  such 
rays  have  been  calculated  and  are  plotted  in  Figure  18.  One  sees 
that  the  regions  A,  B,  and  C  in  Figure  12  are  not  real  shadow  zones 
but  are  interlaced  with  these  narrow  bundles  of  rays.  It  must  be 
observed  that  this  effect  would  not  appear  if  the  speed  of  sound  vs 
depth  profile  had  a  continuous  first  derivative  at  z  =  60  m. 

In  Figure  13  the  deep  sound  channel  with  the  oscillating  rays 
slopes  downwards  with  increasing  range.  The  shadow  zones  A  and  B 
are  substantially  enlarged  as  compared  with  the  corresponding 
shadow  zones  in  Figure  10,  The  range  of  A  has  increased  by  43%. 

Figure  14  shows  that  the  thermal  mixing  zone  has  formed  two 
new  shadow  zones  near  the  surface  —  regions  C  and  D  —  which  are 
not  present  in  Figure  11. 

Figures  15  through  17  show  the  corresponding  ray  paths  with 
the  souhd  source  placed  at  the  point  of  inflection  of  the  thermal 
mixing  zone  (y  =  30  km) .  Here  it  will  be  noticed  that  the  depres¬ 
sion  of  the  rays  is  even  more  pronounced  than  in  Figures  12-14, 
because  now  the  source  is  at  the  position  where  the  effect  of  the 
mixing  zone  is  the  greatest.  Comparison  of  Figures  16  and  10 
shows  an  increase  of  517,  in  the  range  of  A. 

In  connection  with  the  shadow  zones  A,  B,  and  C  in 
Figure  15,  the  same  comments  as  for  Figure  12  apply. 
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Approximate  Equation  for  Total  Angle  of  Deviation  for  a 
Ray  Near  a  Heat  Source 

In  general  the  curvature  of  a  line  is  given  by: 


K  =  d8  - 
ds 


d“=r 

ds* 


From  the  eik.onal  equation  (1)  there  follows: 


d2r’  =  l( Vn  _  dn  dr  \ 
ds2  n\  ds  ds  f 

From  (A-2)  the  curvature  of  a  sound  ray  becomes: 


"4/ws  -  (Vn  •■■#)'  • 


Vn 


or: 


K  =  J  1  -  cos8 1 


where  is  the  angle  between  Vp  and  —  . 

ds 

Equation  (25)  yields: 


Vn  =  -  2kpe"P  (Vp)  , 


(A-2) 


(A-3) 


(A-4) 


where  p  is  the  normalized  radius  p8  =_  §2  +  T]3 

In.  this  problem  the  direction  of  Vn  is  given  by  the  unit  vector 

The  ray  has  the  direction 
The  angle  f  is  given  by: 


Vp  with  direction  cosines  -2  and  -3 

P  P 

M  and  521 


cosines 


ds  ds 
cos  t  = 


l/gdS  + 

p\  'ds  +  'ds 


(A-5) 
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The  curvature  (A-3)  becomes 


2ke“P  /  ' 

- T~nS  /  P"  T 

1-f ke"  P  v 


(A-6) 


The  total  angle  of  deviation  A0  between  the  initial  and  final  direc- 
tions  of  the  ray  is  given  by  the  line  integral 

CQ 

'  A0  =  J  Kds  .  (Ai7) 

s ~Sq 

The  equation  of  the  ray  is  not  known  analytically,  therefore 

this  integral  is  evaluated  with  the  following  approximations,, 

-o2 

Since  k  <  <  1,  ke  K  can  be  neglected  in  the  denominator  of  (A-6). 
Also,  —  <  c  ~~  ,  Therefore  the  approximations  —  0,  T]  T1('~ob)  , 

P3  ^  5s  +  [Tl(“<:n)  l3 ,  —  ~  1  and  ds  *“  d§  can  be  used.  With  these 
approximations  equation  (A- 7)  becomes 

A8  -  2kT|(~eo)e“  d§  , 

—  CO 

which  reduces  to  equation  (29)  in  the  text: 

AG  =  2*/tt  k  Tl(-®)e“t1K"“)^ 
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Typical  Temperature  vs.  Depth  Data  for,  the  Ocean  . 

In  the  table  below  are  listed  the  depths  together  with  the 
corresponding  temperatures  in  Figure  8,  between  which  the  linear 


interpolations  of  h(z)  in  equation  (31)  are  calculated. 

Depth  z  Temp  T 

in  meters  in  C 


0 

60 

200 

250 

300 

350 

500 

600 

800 

1000 

1500 

2000 

2500 

3000 


12.00 

12.00 

7.83 

6.72 

5.89 

5.33 

4.50 
4.06 

3.50 

3.11 
2.67 
2.28 

2.11 
2.00 


5000 


2.00 
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Figure  1:  Ray  paths  near  a  heat  source  with  T0  =  10 °C, 
(*  indicates  position  of  heat  source) 


Temperature  T  in  °C 


1420  1500  1560 

Speed  of  Sound  V  in  m/sec 


Figure  8:  Temperature  T  and  speed  of  sound  V  vs,  depth  z  at  y 


Horizontal  Range  y  in  km. 

Figure  11:  Ray  paths  in  an  ocean  with  the  speed  of  sound  a.  function  of  depth  only 
Sound  source  at  z  =  3000  m. 


Horizontal  Range  y  in  km. 

Figure  12 :  Ray  jpaths  in  ian  ocean  with  a  thermal  mixing  zone .  Sound  source  at  z  *  0  i 
y  «  0  m.  (Arrow  indicates  inflection,  point  of  the  thermal  mixing  zone.; 


o 


Figure  13;  Ray  patte  in  an  ocean  with  a  thermal  mixing  zone.  Sound  source  at 

z  =  600  m,  y  =  0  m.  (Arrow  indicates  inflection  point  of  the  thermal 
mixing  zone . ) 


Figure  16s  Ray  paths  in  an  ocean  with  a  thermal  mixing  zone.  Sound  source  at 

z  **  600  m,  y  =  30  km.  (Arrow  indicates  inflection  point  of  the  thermal 
mixing  zone.) 


Figure  17:  Ray  paths  in  an  ocean  with  a  thermal  mixing  zone.  Sound  source  at 

z  =  3000  m,  y  =  30  km.;  (Arrow  indicates  inflection  point  of  the  thermal 
mixing  zone.) 


ro[  uf  z  ipdag 


Figure  18:  Effect  of  the  thermal  mixing  zone  on  shadow  zones  with  the  sound  source; 

at  the  surface.  Five  rays  are  shown  interlacing  the  shadow  zones  A,  B, 
and  C  of  Figure  12. 
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